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1 Introduction 

The three dimensional Navier-Stokes equations are a not yet completely under- 
stood mathematical problem, in the sense that there is no proof of uniqueness 
of solutions in the spaces where existence is proved. This mathematical prob- 
lem has been investigated since long, also in connection with the analysis of how 
good are the Navier-Stokes equations to model turbulence. Some attempts have 
been made to overcome the problem of uniqueness, introducing some modifica- 
tion in the Navier-Stokes equations. In this paper we are concerned with the 
model proposed by Prouse in [TS]. Here, we study a stochastic version of this 
problem, as explained below. As soon as a stochastic equation is introduced, 
statistical properties typical of turbulence can be investigated. 

We remind that nowadays there are many results on stochastic three dimen- 
sional Navier-Stokes equations (see, among the others, [T], [2], [3], @], [5], [B], 
[3], [5], [H], [IH]); however, the uniqueness problem is not solved also in the 
stochastic framework. 
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Let us consider the partial differential equations of Navier-Stokes type 



du + [-A$(u) + (it • V) u + Vp - Vdiv $(u)] dt = G (u) dw 



div u = 
u\ t =o = u a 



(1) 



where, for t > and x G T C R 3 , u — u{t,x) is the velocity vector field, 
p = p(t,x) the pressure field; v > the viscosity coefficient. G is an operator 
acting on the noise and on the velocity; the vector function $ : R 3 — > R 3 is 
defined as follows: 



where a-i > a\ > and b > 4. 

$ describes the nonlinear relationship between the stress tensor and the defor- 
mation velocity tensor, as explained in [15]. When this relationship is linear, 
then = vu and (JTJ) are the usual Navier-Stokes equations for an homoge- 
neous incompressible viscous fluid with random forcing term; indeed, the first 
equation becomes 



For problem lfl ]) -([2 )) . in Section [3] we prove a result on existence and unique- 
ness of martingale solutions (Theorem [6|) and on existence of stationary mar- 
tingale solutions (Theorem [9|) . 

In Section 21 another model with $(it) = ;/|u| 4 u in fl} is investigated; anal- 
ysis on existence of martingale solutions and stationary martingale solutions is 
presented (Theorem 1 10[) . Moreover, introducing a scaling transformation sug- 
gested by turbulence theory, some insights in the behaviour of the function 
structure of any order p are shown (Claim [T4|. 

Preliminaries are in Section [21 auxiliary results are in the two Appendixes. 

2 Notations and preliminaries 

Let the spatial domain be a torus, i.e. the spatial variable x belongs to T = 
[0,i] 3 and periodic boundary conditions are assumed. L 2 is defined as the 
space of vector fields u : T — > R 3 with L 2 (T)-components. For every a > and 
p > 1, W a ' p is the space of fields u £ L p with components in the Sobolev space 
W a > p (T). For a < 0, W Q ^ is the dual space of W- Q < p ' with 1 + ^=1. Set 
H"=W Q - 2 . 

We introduce the classical spaces for the Navier-Stokes equations (see, e.g., 
[T?]). T)°° is defined as the space of infinitely differentiable divergence free 
periodic fields u on T, with zero mean (J T u(x)dx = 0). 

Let H be the closure of 2?°° in the L 2 -topology; it is the space of all fields 



$(w) = er(|u|)w with 



a € ^([O.oo)) 

<r(0 > v > 0, o'(i) > 

ai^' 1 < cr(0 < a 2 £ fc ~ 1 when £ > K 



(2) 



du + [—vAu + (u-X7) u + X7p] dt = G (u) dw 
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ueL 2 such that divu = 0, u ■ n on the boundary is periodic, j T u (x) dx = 0. 
We endow H with the inner product 

{u 1 v) H = -^ j u(x) ■ v(x)dx 

and the associated norm \-\ H . 

Let V (resp. D(A)) be the closure of X>°° in the H^topology (resp. H 2 - 
topology) ; it is the space of divergence free, zero mean, periodic elements of H 1 
(resp. of H 2 ). The spaces V and D{A) are dense and compactly embedded in 
H (Rellich theorem). Due to the zero mean condition we also have 

j \Du(x)\ 2 dx> A J \u(x)\ 2 dx 

for every u e V, for some positive constant A (Poincare inequality). Here 

\Du(x)\ 2 =Y. 
inner product 



\Du(x)\ z — J2i j=i (dj u i( x )) 2 ( an( i dj — jUt)- So we ma y endow V with the 



(u,v) v — \_] / 9jUi(x) djVi(x) dx 



and the associated norm ||-|| v . 

Let A : D{A) C H — > H be the operator = — Ait (componentwise). 
There is a complete orthonormal system in H made by the eigenvectors h k .j of 
the operator A (Ah^j = \k.jh k ,j)- Since the spatial domain is the torus, we 
know the expressions of these eigenvectors with their eigenvalues. Indeed, let 
k = (ki,k2,ks) with integer components, i.e. fceZ 3 . We denote by 1? + the half 
space of Z 3 defined as = {k x > 0} U {k x = 0, k 2 > 0} U {k x = 0, k 2 = 0, k 3 > 0}. 
Then for any k £ Z+, there exist two unit vectors Vk,i and v k ,2, orthogonal to 
each other and belonging to the plane orthogonal to k. Then the (four sequences 
of) eigenvectors are 

hk,i(x) = jj^Vks cos(^fc • x) , h ka (x) = j^v k ,2 cos(^fc • x) 

hkA x ) = JJj2 V k,i sm{^-k ■ x) , h kA (x) = jT^v k , 2 sin(^-k ■ x) 
with eigenvalues 

Afe,i = A fci2 = Afc,3 = XkA = ^ \k\ 2 

for any k G Zj_. 

Hence, i7 = span{/ifc ; j : j = 1, 2, 3, 4 and fc G Z 3 } and we set i?„ = span{h k ,j ■ 
j = 1,2,3,4 and k G Z 3 ^, |fc| < n}; moreover, we denote by 7r„ the projection 
operator from H (or any subspace, as V or D(A)) onto H n . The operators A 
and 7r„ commute. 
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We may take the Poincare constant A above equal to (2ir) 2 / L 2 (the first eigen- 
value of A). Notice that we have 

(Au,u) H = \\u\\ 2 v 



for every u £ D(A), so in particular 



(Au, u) H > 



M 2 

L 2 



u 



2 

H ■ 



Let V be the dual of V with respect to the iJ-norm; with proper identifica- 
tions we have V C H C V with continuous injections, and the scalar product 
(■,-) H extends to the dual pairing (•, -) vv , between V and V' and to the dual 

pairing (•, -) L9 L ,/ between L 9 and h q (- + = 1). 

Let B (•, •) : V X V — ► V be the bilinear operator defined as 

(w,5(tt,'ij))y F = J Ui(diVj)wjdx (3) 

for every u,v,w € V. By the incompressibility condition, we have 

(B («,«),«) =0, (B(u, «),«;) =-<fl («,«;),«) (4) 
Using the latter relationship, by Holder inequality we estimate 

\B(u,u)\v= sup |(fl(«,tt),^)|<H£ 4 (5) 

IIV'II V'<1 

We list here a number of inequalities. 
Lemma 1 

(A<f>(u),u) H > HMIv 

Proof. We have 

3 

= (<f>(u),u) v = / V" [dk$i{u)]dkUidx 
Jt 'i, fc =i 

The estimate on ^ i fc=1 [5fe$i(u)] SfcUi comes from [15]. ■ 
Lemma 2 

($(«(!)) - $( W ( 2 )), U (!) - U ^) H > V\U^ - U^\ 2 H 

Proof. The proof is by [15]. We rewrite it here, because we shall need it in 
Section [4] 
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Set cr(M) =u + a(\u\) with a' > 0. Then 



r(\uW\)t 



(1) 



r(\u^\)r 



(2)1 



,(1) 



>a(\u^\)\u^\ 2 + a(\u^\)\u^\ 2 

- ki^Di^i 2 - ~*(l« (2) l)k 2) l 2 - kl« (1) l)k 2) l 2 - kk 2) l)l« (1) l 2 



2 
1 

~ 2 
> 

Hence 



ia(|««|)| U W| 2 + la(\ u W\)\uW\ 2 - la(\uM\)\ U W\ 2 - ^(| U ( 2 )|)|^)| 



2 u 1/1 1 2 
-a(|u (2) |)l - |u (2) | 



,(2)1 



= u\u^-u^\ 
>v\u^-u^\ 



\a(\u {1) \)u (1) 



a(\u^\)u^\ 



,(2)1 



Next lemma is crucial to prove uniqueness. Notice that the regularity u G 
L 5 (0,T;L 5 ) is needed here. The weak solutions of the Navier-Stokes equations 
(deterministic or stochastic), which are known to exist, are not proved to have 
such a regularity; here the modified term $ (with b > 4) plays its role. We 
remind that Prodi [13] proved uniqueness for the deterministic three dimensional 

1q 

Navier-Stokes equations, if u G Li- 3 (0,T;L 9 ) for some 3 < q < oo. For q = 5 
the required regularity is u G L 5 (0,T;L 5 ) and this implies uniqueness also in 
the Prouse model (see [15]). 

Lemma 3 If u G L 5 and v G H , then for any v > 

\(B(u,v),A- 1 n m v)\ < ^\v\ 2 H + C B \u\t 5 \-K m v\ 2 v , 

{{Bty^A-^mV)] < ^\v\ 2 H + C B \u\l 5 \TT m v\ 2 v , 
for some positive constant C b ■ 

Proof. In [15] . there is a very similar lemma, but with v instead of n m v (here 
we consider any finite projection operator n m ). Following the lines of that proof, 
we get our result. ■ 

Properties of G 

Let G : H —> L (H) be a mapping with the properties 



\\G(u)\\hs(h) < AoHlr+P 



(6) 
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and 

A-^[G(v) - G{z)\ ' < L G \v z\l, (7) 

Ho\H) 

Here H^H^s^) is the Hilbert-Schmidt norm of an operator in H, defined as 



Now, we project equation |T]) onto the space of divergence free vectors fields; 
both the V-terms desappear, as when we deal with the Navier-Stokes equations 
(see, e.g., [IB]). Then, we obtain an evolution equation (still formally), which 
with our notations is 

du + [A$(u) + B( u, u)] dt — G (u) dw, u(0) = uq (8) 

From now on, <i> will be assumed to satisfy Q for a given b > 4. 

The rigorous interpretation of this equation will be given in the sequel, but for 

the time being let us at least write it in weak form 

(ut,ip) H + J mu s ),Aifj) hl+ i^ 1+b ds-J^ (B(u s ,tp) ,u s ) L | L4 ds 

= {uo,^) H + (G( u s)dw s ,ip) H (9) 



with tfj e V°° and < t < oo. 

We assume that w is a cylindrical Wiener process in H (see, e.g., [7])- We can 
represent it as follows. Suppose we are given a Brownian stochastic basis, i.e. 
a probability space (W,JT, Q), a filtration (Tt)t>o and a sequence {/3k,j(t)}k,j 
of independent Brownian motions on (w, J 7 , (J r t) t>0 , Q^j . Namely, for k G 1? + 
and j = 1,2,3,4, the real valued processes f3k,j(t) are independent, adapted to 
(•^*)t>0' continuous for t > and null at t — 0, with increments 0kj (t) — 0kj (s) 
that are N (0,t — s)-distributed and independent of T s . Then 

4 

j=i ka%\ 

is a cylindrical Wiener process in H. 

The convergence of this series requires proper distributional topologies. The 
stochastic integral in equation © is well defined under the Hilbert-Schmidt 
assumption made on G (see [7] for details). 
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3 Well posedness 



3.1 Concepts of solution 

Consider the abstract (formal) stochastic evolution equation (JSj) and its weak 
formulation over test functions (O . We have 

$(u s ),^V) L i + i L1+b |rfs < \$(u s \ 1+ t \A^\ u+b ds 

< l (1 + \u s \ u+b )ds 
Jo 

because 

l*( u )£T+\ = J \<5>(u(x))\ 1+ h {lu{x)l < K} dx + mu(x))\ 1+ h {lu(x)l>K} dx 

by @ 1 f i f , 

< K 1+ ^ J^\a(\u(x)\)\ 1+ h {lu(x)l < K} dx + Jja 2 \u(x)\ b ) lH dx 

<C$(1 + J \u(x)\ 1+b dx) 

(11) 

since a € C 1 implies that a is bounded on [0, K\. Then, in equation (J9j> the term 
Jo ($(u s ),Aip)ds is well defined for functions u that live in L 1+b (Q, T; L 1+h ). 
T > 0. 
Moreover, 



(B(u s ,ip),u s ) ]L i^ 4 ds<J^ \u s \ hi \\ip\\ v ds < J \u s \ Li ds (12) 

Hence, in equation @ the term J^(B(u s> , ip),u s }ds is well defined for functions 
u that live in L 2 (0,T;L 4 ). 

We conclude, in both cases, that given b > 4 the regularity u £ L 1+b (0 7 T; L 1+b ) 
is enough to define these quantities. Moreover, from now on the duality pairing 
for these two terms has to be understood in the sense above specified (as written 
also in equation ((9])). 

As in the deterministic case, strong continuity of trajectories in H is an open 
problem. There will be strong continuity in weaker spaces (like W -2-9 ' 1 " 1 "*), 
and a uniform bound in H. Let H a be the space H with the weak topology. 
Since 

C([0,T];W- 2 - e ' 1+ ^) n L°°(0,T; H) C C([0, T]; H a ) 

then the trajectories of the solutions will be at least weakly continuous in H 
(see [IS] pg. 263). 

Given a separable Banach space W (it will be W = W~ 2 ~ e,1+ t ), let us set 

O = C([0,oo);W') 



7 



and denote by (£t) t>0 the canonical process (£ t (lo) = ui t ), by F the Borel a- 
algebra in f2 and by F t the a-algebra generated by the events (£ s <E A) with 
se [0,i] andieB(W'). 

Definition 4 (solution to the martingale problem) Given a probability mea- 
sure Ho on H , we say that a probability measure P on (Q, F) is a solution of 
the martingale problem associated to equation Q) with initial law /io if 

[MP1] for every T > 

p[ sup 161*+ [ T \\a 2 v ds+ f T \^lda< oo) =1 
\te[o,T] Jo Jo J 

[MP2J for every -0 £ T>°° the process Mf defined P-a.s on (fi, F) as 



JO JO 

is square integrable and (^hlf ,F t ,P^j is a continuous martingale with 
quadratic variation 

M t = [ \G(^\ 2 H ds 
Jo 

[MPS] /io = IIoP, where ITo denotes the restriction on Fq. 

Remark 5 j4 solution of the martingale problem is also a weak solution. The 
definition of weak solution is as follows: there exists a Brownian stochastic basis 

t>0;k.j ) an d a -valued process u on (W, T . Q) such 

that 

[WM1] u is a continuous adapted process in W and 

u{.,u) eL°° (0, T; H) n L 2 (0, T; V") n L 1+fc (0, T; L 1+b ) Q-o.a. 

/or every T > 

/WM2/ (0j is satisfied Q-a.s. 

[WM3] u(0) has law fj, . 

Finally, in this context, we call strong solution a process u satisfying the 
three above properties on any a priori given stochastic basis. 
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3.2 Main result 

Theorem 6 Let p be a measure on H such that m p := J H \v\jj /j,(dv) < oo for 
some p > 2. Then there exists one and only one solution to the martingale 
problem @ with initial condition \x. 

Moreover, two strong solutions on the same Brownian stochastic basis coin- 
cide a.s. 

Proof. Step 1 (Galerkin approximations). Let 

(Ft)t>0,Q,(PkAt))t>0;k tj ) 

be a Brownian stochastic basis supporting also an JFo-measurable r.v. uq : W — ► 
H with law p. For every n, let Uq := 7r„uo and consider the Galerkin system 

du? + L4$«) + 7r„B«, u*)] dt = ir n G «) dw t , u n (Q) = v% (13) 

obtained by applying the projection operator tt„ to both sides of equation (JSJ). 
(Notice that ir n A$(u?) = A$(u?).) 

Equation (|13p is a stochastic ordinary equation in the finite-dimensional Hilbert 
space H n . 

Local existence and uniqueness (on a random time interval) is classical, since the 
nonlinearities are locally Lipschitz continuous (see, e.g., [H]). Global existence 
is then a consequence of the a priori estimates given in Appendix 1. There, 
defined Tjj = inf {t > : \u\ L \ 2 H = R) we shall prove that, for any T > 

E sup [uSv^l* < Ci (14) 

0<t<T H 

E [ T \K AT J 2 v l {s<TR] ds < C 2 (15) 
Jo 

E f T \u^ AT „\lt^ {s< ^}ds<C 3 (16) 

** 

for some positive constants G\ — Ci(p,T, Ao, p,m p ), C% — Ci(T, Ao, p, TO2), 
C3 = 03(0.1, G\, C2), independent of n and R. 

Now, assume first that the initial velocity is bounded: \uq\h < K. Take 
R > K; so > Q-a.s.. The solution u™ to the Galerkin system (|13|) is 
defined at least in the time interval [0, rjj). Since we know from (fT4|l that 

£ sup |< AT „| 2 ff <Ca 
te[o,T] 

for some constant C\ — C^ p independent of n and R, we have 
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for T > fixed. Moreover 

Q(r£ <T) = E1 {T , <T} = ±E (i {tR< t}K AtR \ 2 h ) 
because \u^ hT ^\ 2 H = R on the set {t^ < T}. Hence 

Notice that t~ > rj| for R> R. Therefore, setting = sup the process u™ 

R>if 

is defined for t £ [0, t^). But we have 

Q(r£ < T) < Q(t£ < T) < ^ Vi? 

Hence 

Q(r* < T) = 

and finally we conclude that w" is a solution for t G [0, T). Since T has been 
chosen arbitrarily, we conclude that the Galerkin solution is defined on any finite 
time interval. 

For a general initial velocity satisfying the assumption of Theorem [6l we 
proceed as follows. Let Wk G T be defined as Wir = {|uo|lf < K}; we have 
Q(UifWif) = Q(|uo|h < oo) = 1. Define uo K as uo on Wjf and otherwise. 
Let it™^- be the unique solution to the Galerkin system (Vt > 0) with initial 
condition uq k . EK>K, then 

Q {W K n = u n tK Vt > 0}} = Q{W K } 

We may uniquely define a process u™^ on W = Ua-Wk as it™^ = u™ x on 
Wjf. Looking at the Galerkin equation in the integral form, it is clear that 
u too so l ves the equation on W' . But Q(W) = 1. Thus we have proved the 
existence of a global solution to the Galerkin system for any initial velocity 
with m p < oo for some p > 2. This solution is a continuous adapted Markov 
process in H n (uniqueness holds for the Galerkin problem; it can be checked 
direclty or obtained as a byproduct of next Step 5). 
Hence we have proved that, for any T < oo 

E sup \v%\ p H < d (HI) 

0<t<T 

E C \\u n s \\ 2 v d S <C 2 (SI) 
Jo 

e [ T K\lt+*ds < C 3 (Hi) 
Jo 
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From these estimates, we also get the following one. Given ip 6 T>°° and 
e e (0, 2), we have 



E 



(-K n G(u")dw s , V) 



2+e 



< E 



4\ 



(ir n G(v%)dw s ,ip) 
= \C\E ( \K n G(u n s )^\ 2 H ds 



by Gaussianity 



\HS{H) 



ds 



2\ 4 



<ch^ / (K|£* + i)ck 

Jo 

Here (and in the following) C denotes different positive constants, independent 
of n. Taking 2 + e < p and bearing in mind (|14|) . we conclude that for any finite 
t 



supE 1 



(ir n G(v%)dw s ,il)) 



2+e 



< oo 



(17) 



Here the limitation e < 2 can be easily removed, but in the sequel it will be 
enough to consider a positive quantity s as small as we want. 

Step 2 (time regularity and reformulation in path space). In view of the 
time regularity, equation (|13| has the form 



where 



It 



A$ «) ds 



Jf = - / 7r n fl(u;X)ds 



X? = / Tt n G(u n a )dw 8 



For the first term we have 



tr 



(o,T;W _2 ' 1 + £) 



Jo L ¥ 

<C(T + [ T \uXt^ds) 
Jo 
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according to (fTTj) . 

For J", using ([5]) we have 

ll^ll^(0,T;V0 < C f \B{u n s ,U n s )\ 2 y,ds 
J 

<c f T w:\ids 

Jo 

<C(T+ f K|[+Us) 
Jo 

Finally, for every q > 1, a G (0, |), T > 0, we have (see, e.g., [5]) 

^ii^"ii^, 9(0 ,t ; «) f hnG{u:)r HS{H) d S 

J 

and by ([5]) and the mean estimates of the previous step we conclude that 

E \\^ n \\w a '"(0,T;H) — C 

(C independent of n and p > 2 as stated in Theorem [SJ) 
Therefore, for a £ (0, |) 

u n e T . W -3,l+i^ + ( Qj J,, y/j + ^a,p ^ T . H ^ 

in mean. Notice that H C F' C W~ 2 ' 1+ i M /1,2 (0, T) C W^ 1+ i(0,T) C 
W^«.i+i(0,T) and PF Q < p (0,T) C W^ 1+ i (0, T). 
We conclude that, in mean 

Under the embedding C -ff , we have that (w")t>o is a continuous adapted 
process in H, so it defines a measure P n on C ([0, ooj; and thus on (f2, F). 
Actually, P n is concentrated on C([0, oo);H n ). For every a G (0, |), T > 0, the 
above estimates may be rewritten as 




and 



E 1 



II" 



< C 5 (T, X ,p,m p ,b) 



(19) 



■ 1+ i(o,T;W~ 2,1+ i; 
for any n. 

Relationships (|14l V(|16l) may be rewritten in a similar way. 
Step 3 (tightness). Use now Chebyshev inequality and (|15l ). ([TO)) . Then, 
given ag (0, |), T > 0, for every e > there is a bounded set £? e such that 



B £ C L 2 (0, T; U) n TU Q ' 1+ > (0, T; 
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and 



infP„(P e ) > l- £ 



The space L 2 (0, T; V)nW a ' 1+ i; (0, T; W^ 2 *^) is compactly embedded in L 1+ i (0, T; 
(see, e.g., Theorem 2.1 in |9]). Hence, for every e > there is a compact set K £ 
such that 

K e C L 1+ ^ (0,T;H) and inf P n (if £ ) > 1 - e 

n 

Now, take any separable Banach space W such that W~ 2 ' 1+ ^ is compactly 
embedded in W; e.g. W = W^ 2 " 9 * 1 ^ for some 9 > 0. Notice that all the 
spaces VF M+ £(0,!T), W l ' 2 {Q,T), W a ' p (0,T) (for ap > 1) are continuously em- 
bedded into C([0,T]). Hence, the space of vectors with the regularity specified 
by Ull) is compactly embedded in C([0,T]; W) (see, e.g., Theorem 2.2 in [9]). 
From the boundedness in the mean of I n in H /1,1+ 's(0,T; W^ 2 ' 1+ ^ ), of J n in 
W 1 ' 2 (0, T; V) and of the law of the Wiener process in W a p (0, T; H) for every 

a € fj> 5^ i again by Chebyshev inequality and compact embedding we obtain 
that for every e > there exists a compact set K' e such that 



K>CC([0,T}. 



and infP n (^)>l-e 



Therefore the family of measures {P n } is tight in (Q,T;H) and in 

C ([0, T] ; W), with their Borel cr-fields. Hence there exists a probability measure 
P on 

C([0,T];W')nl 1+ 5 (0,T;H) 

that is the weak limit in such spaces of a subsequence {Pn k }- 

Step 4 (P is a solution to the martingale problem). From the uniform 
estimates on {P n J in i 2 (0,T;^), L°° (0,T;H) and L 1+b (0, T; L 1+b ) we may 
deduce that P gives probability one to each one of these spaces and has bounds 
in the mean similar to those uniform of P„ fc . This way we have checked property 
[MP1] in the definition of solution to the martingale problem. 

Concerning [MP3], we have P nk — > P as weak convergence of probability 
measures on C([0,T]; W); in particular noP„ fc — > n P as probability measures 
on W. But noP„ fc is the law of Tr nk uo, which converges to /i since n rik uo 
converges Q-a.s. to uo- Hence HoP is fx. 

Finally, let us check property [MP2]. We proceed as in [7] (Sec. 8.4) or in 

m- 

Given tp £ T>°° , we have to prove that for every t > s > and every bounded 
P s -measurable random variable Z, we have 



E 1 



E P [(Mt) 2 ' 
(Mf-Mf)z 
[{Mff-, S ])Z 



< oo 

= 
= 
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where ft : = J |G(^ a )^|^ds. Defined 



for the measure P„ fc we know (see, e.g., [7] Sec 8.4 ) that (^Mf ,nk ,F t ,P n 
square integrable martingale with quadratic variation 



[M^' nk ] t = $ k = 



\^n k G^ s )ip\ 2 H ds 



Thus 



([{Mp nk ) 



,\2 _ ^rifcj 



Z 



(M?' nk - Mf' 
[(Af^' n *) 2 -«?*]) Z 



= 
= 



Moreover, by (fTT|) we know that there exists some e > such that 

2+e 



sup£" 

k 



M i>,n h 



< oo 



(20) 
(21) 

(22) 



Now, let us consider the limit as k — > oo. 
We know that P nk converges weakly to P; then by Skorohod theorem there exists 
a stochastic basis (Q, F, F t , P) and, on this basis, there exist L 1+ £ (0,T; H) n 
C([0,T];W) -valued random variables u,u nk such that u has the same law of 



has the same law of 



and 



C([0,T];W')-norm. 
Define 



P-a.s. in the L 1+ i (0,T; H) H 



M, 



{u n t \i>)H-{u n \i>)H- Ai/>)da 



Then ([20| - ([22|) hold true (with the obvious change of notation). 
If we prove that Mf' nh — > Mf P-a.s. as k — » oo, then by the equiboundedness 
relationship ((22]) we obtain that M"/' n * -> M/ in L^O, P) and in L 2 (f), P) and 
f" fc — * ft in L 1 (0, P). This concludes our proof. So, we have to prove a P-a.s. 
convergence for each term in the definition of Mf' Uk . 
It is trivial that P-a.s. 

(Ut",1p)H ~ (uo k ,ip)H -> (u t) ip)H ~ (uo,iP)h 
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Notice that there appears the scalar product in H and not the duality pairing 
(lit, V')w',W) because the limit process u belongs to C([0, T];H a ) with probability 
one. 

Moreover, there exists a subsequence (we do not write that we consider a 
subsequence, since we shall pass through subsequences a few times from now 
on) such that 



P-a.s. u" h (x) 
We also have, for any k 

E 



,(x) for a.e. (s,x) G [0,T] x T 



?,rik 1 1+6 



■it 



a \l. 1 + b 



ds < C 3 



Keeping in mind pip and ([5]), it follows that $(u™ fc ) is equibounded inL 1+ f>(J7x 
[0,T] x T) and (B(u rik , 7r„ fc ip), u nk ) is equibounded in L j t ± (0 x [0,T]) respec- 
tively. First, we get that $(it™' c (a;)) — > $(m s (s)) P-a.s. and for a.e. (s, a;). By 
the equiboundedness of $(-u" fc (a;)) in L 1+ i(£lx [0,T] xT) it follows that <&{u nk ) 
converges to in x [0,T] x T); we get 

E [ (<5>{u™ k ),A4>)ds -> £ / 
Jo Jo 

Hence a subsequence of / *($(u™ fc ), ^4V)^ S converges P-a.s. 

On the other hand, another (sub)subsequence can be extracted so that 

P — a.s. u™ k —> u s in L 4 for a.e. s 

Then, by triangle inequality (B(u™ k , 7r„ fc t/>), u 7 s lk ) — > (B(u s ,ip),u s ) P-a.s. and 
for a.e. s. By the equiboundedness of (B(u" k , 7r„ fc -0) , w" fc ) in L^Jl x [0,T]), 
we conclude as above that there exists a subsequence of f*(B(u™ k ,?p),u™ k )ds 
converging P-a.s. . 

Considering the convergence of a suitable subsequence (the last extracted), 
we get that (|20l) - (p?!?|) in the limit allow to conclude the proof. 

Step 5 (uniqueness). Let i = 1,2, be two strong solutions on the same 
Brownian stochastic basis. We are going to prove pathwise uniqueness, which 
implies uniqueness of martingale solutions. 

Let 



6t = 2Cb 



4 2) , 

,(1)|5 



,(2)|5 



with Cb as in Lemma [3] and Lq as in 0. 

We have 



de-SS'-*>\v?\* 



he 
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By Ito formula, the last differential is 



d\v r ' 



2(7r m [$( u «) - $(« (2) )], v m ) H dt - 2(ir m [B(u^, v) + B(v, u^)), A~ l v m ) H dt 



+ 2(n m [G(u^) - G(u< 2 >)] dw, A~ 1 v m ) H + 7: m A^ 2 [G(u^) - G(u^)} 
We estimate some terms as follows. By Lemma [3] 

2\(i: m [B(u {1 \v) +B(y,uW)],A- 1 v m )\ = 2\{B(u^, v) + B(v, A^v™)] 

<p\v\ 2 H + 2C B {\ui 1) \l + \u^\l 

By Lemma [5] 

2(n m [t>(u^) - ${uW)],v m ) H = 2{<S>{v m u^) - <5>{Tt m u^),v m ) H - 2e 



dt 



HS(H) 



> 2v\v m \ 2 H -2e r ' 



li 5 (0,T;L 5 ) 



where/ T |erl^<C(| 
By © 

n m A-^[G(u^)-G(u^)} 



/.( 2 )| 



li 5 (0 



5 ,) and lim / T ef dt = 0. 



< 



HS(H) 



A- 1 /2[ G(u (l))_ G(u (2) )] 



HS(H) 



< L G \v\y, 



Now, we integrate in time equation (|23p and use the above estimates, ob- 
taining 



X e ° ds \v%\ 2 v ,+v / e-/o t ».*(2|« t "*|ff-|« t |ff)d« 



< 



K\ 2 v ,+2 e-ti e ° d °e?dt+ L G e-ti e ° ds (\v t \ 2 vl -\v?\ 2 v ,)dt 



+ 2 I (n m {G(uP)-G(uY')]dw t ,A- l vr) H 
Jo 

We can take the limit as m — > oo in every term. We get 



(2)> 







< h| 2 v ,+2 / ([G(^ 1) )-G(^ 2, )]du; t ,^- 1 « t ) ff 



(2h 



Hence 
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When the initial conditions of coincide, we deduce 



- So Ssds 



\v t \ 2 H dt = 



with probability one. Since f Q 9 s ds < oo a.s., we have v = consider- 
ing v as a measurable function of t with values in H. This implies that with 
probability one w 1 ' = vS 2 \ where the equality holds in L°°(0, T; H). ■ 

3.3 Markov and Feller property 

Lemma 7 Let Uq,Uq be initial data satisfying the assumption of Theorem & 
and let (u™)t>o ond (u t )t>o be the corresponding strong solutions on the same 
given Brownian stochastic basis. 

If E\uq — uo\y, — > 0, then for every T > 0, (it™)*>o converges to (ut)t>o in prob- 
ability on [0,T] x £1 in the topology of H, and converges to ut in probability 
on in the topology of V' . 

Proof. We proceed as in the previous Step 5 to get the following estimate 



E 



ut\ v > 











+ vE 


/ e-J 
Jo 


° e > ds \u?-u t \ 2 H dt 


<E\u%-u \ 2 v 



where 



= 2C B [\u?\l + \u t \l] +L G 



Since 9 s ds < oo with probability one, we get the result. ■ 

Theorem 8 The strong solutions of equation (0) on a given Brownian stochas- 
tic basis define a Markov process in H with the Feller property in V . 

Proof. Denote by u(t; y) the solution at time t which started at time from y. 

Given t > the dynamics y i— > u(t; y) is uniquely defined in H; hence the 
Markov property is inherited by u from the Galerkin approximations u n . 

The process solution enjoys the Feller property if 

Eg(u(t; z)) Eg(u(t; y)) as z — > y in V' 

for any t > 0,g € Cb(V). For this it is enough the convergence in probability: 
u(t; z) — > u(t; y) as z — > y in V' . But, as in Lemma [7] (now the initial data are 
deterministic), we know that 



E 



e-ti e ° ds \u{t;z)-u{t;y)\ v , 



Then, we conclude as before that \u(t; z) 
\z-y\v -> 0. ■ 



< \z 
u(t;y)\v< 



y\ 2 v> 



in probability as 
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3.4 Stationary solutions 

As in [9], existence of stationary solutions is obtained in the limit, showing first 
that the Galerkin problem has at least one stationary solution. Our result is 
the following 

Theorem 9 Assume that , { > Aq. Then equation (jHJ) has a stationary 
solution. 

Proof. Let us consider 

du™ + + 7r„B(u™, u™)] alt = w n G «) dw t , = 

We use estimates from Appendix 1. By (f29|) , using > x"M# we 

±E\u^ H+P ^E\uir H < |p(p-l) [A S|<|^ + pi?|<|^ 2 

< ip(p-l)(Ao + e)i?K^ + C(e,p,p) 

for some positive e. 

If 2z/ v J > Ao, then there exist p > 2 and £ > such that pv^j4~ > ^p(p ~ 
l)(Ao + e). Therefore there exists a > such that 

^£K |* + a£K Iff < P) with «o = 0; 

by Gronwall Lemma we get 

SK|^<C 6 Vi>0,Vn>l 

Hence, the family of random variables {u™}t>o is tight in H n , Notice that the 
Galerkin problem is Feller in H n . Then, by the Krylov-Bogoliubov method we 
get that there exists a stationary solution (whose law we denote by /x„) for the 
Galerkin equation. 

Now, consider the Galerkin problem with initial velocity of law /i n and de- 
notes the law of the solution by P n (a probability measure on C([0, oo); W). 
We have 

E P "\i Q \ p H <C & Vn>l 
The corresponding solution P n is a stationary process in H n , i.e. 



P n Ut+ f A$(e,)ds+ f 7r n B(e„C a )ds = Cr+ j 7r„Gfe 

I. J r <J v J T 



)dw ( 



and 

Pn(6) = Pndr) 

for any < r < t < oo. 
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Now we proceed as in [5]. Endow L 1+ & (0, oo; H) with the distance 
(u, v) = —r (\u — v\ ,,i A ll 

fe=l 



and C([0, oo];W) with the distance 

oo ^ 

doo(w,«) = J! U ~~ w lo([Q,*];W) A l) 
fe=l 

The convergence with respect to d 1+ i + <ioo is equivalent with the convergence 
on every finite time interval. We come back to the bounds (|14l ). (|15l ). I|16l ) and 
(fl"§|) , to notice that they hold true because they depend only on E\uq\ p . Thus 
we get tightness on every finite interval; we pass to the limit for a subsequence 
and get the limit process P which is stationary, since the P n are so. It can be 
shown as before that P is a martingale solution to equation jS]). 

Defined the Markov semigroup P* acting on the space of Borel bounded 
functions Bb(H) as Pt4>(y) = E(j>{u(t\ y)), we get that the law /j, of this station- 
ary solution is an invariant measure, in the sense that J Pj^cf/i = J (pdfj, for 
any G B b {H) and t > 0. ■ 



4 The case <f>(u) = v\u\^u 

Instead of ©, let us assume that 

= i/\u\ 4 u 

This corresponds to the case 6 = 5 with the nonlinearity acting everywhere. 
The interest in this model will be explained in Subsection 4.2. 

We can analyze this model as done in the previous section, with few changes. 
Mainly, the solution will live in L°°(0, T; H) n£ 6 (0, T; X), where X is the closure 
of T>°° w.r.t. the norm 

Mx = ^jf {|u(aO|*|Vu(aO| 2 + 4| u (x)| 2 • d % u{x)] 2 }dx^j 

Notice that the term J Q ($(u s ), Atp) H ds in the equation is well defined, since 
\{<f>(u a ),Aip)\ds < \Aip\w f |$(u s )| L6/5 ds = \Aip\ h 6V [ \u s \l 6 ds 



The last integral is well defined for functions u £ L 5 (0,T;L 6 ). But, if u € X, 
then it G L 6 by Theorem [T5l in Appendix 2. 
We have the following result 
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Theorem 10 Let = ^|w| 4 u. 

Let [i be a measure on H such that m p := J H \v\ p H [i{dv) < oo for some p > 2. 
XTien there exists at least one solution to the martingale problem © with initial 
condition fi, assuming that condition [MP1] in Definition 5 is replaced with 



P\ sup \t t \ H + f 
\te[o,T] Jo 



£ s \\lds < oo = 1 



Moreover, if IvCx > Ao (with \u\^ > Cx\u\ 6 H ), then there exists a stationary 
solution. 

Proof. Let us check step by step how our previous proof (for 6=5) can be 

adapted to handle this model. 

Step 1 Instead of Lemma (QJ , we use 

{A${u),v) H = v\u\i 

Hence if we apply Ito formula (for p > 2) to |u"at" \ v Hl we get (JUJ); however (|15[) 
is replaced by 

i-T 



Ej \K AT jihs<r R }d S <C' 2 

Jo 



p^|) is a consequence of the latter relationship, since XcL 6 . 
Step 2 The estimates are still valid: 

sup^Kl^i (0)T;W - 2 , f) < oo for < a < \ 

and 

sup£|M n | L 6 (0iT . X ) < oo 

n 

Step 3 What we need is a compact embedding, which is given in Theorem 
QH Thus, by Theorem 2.1 in 9 the space L 6 (0,T;X) n W a >s(0,T; W" 2 -!) is 
compactly embedded in L^(0, T;L 6 ). Therefore the family of measures {P n } 
is tight in Li(0,T;JL 6 ) and in C([0,T];W'), chosen W such that W~ 2 'f is 
compactly embedded in W'. 

Step 4 The remaining part of the proof for the existence holds true. 

Step 5 As far as the uniqueness is concerned, Lemma[2]has to be replaced with 

- $(w (2) ),u (1) - u^) H > 

This comes from its proof, when we put a = a. The above inequality is not 
enough to get uniqueness. Other estimates failed to be useful so far and unique- 
ness is an open problem. 

Stationary martingale solutions. We consider the sequence of Galerkin solu- 
tions {it™} n >i, all with zero initial velocity. From the estimates in the Appendix 
1, we get 

jE\u^ H +pvE[\u^- 2 \u^\l] < lp(p-l)\ E\u?\P + \p{p-\) P E\u^- 2 
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By the embeddings X C L 6 n H, we get > C x \ u\^j\ thus 

j t EK\% + P vC x E{\uIT h 2 K\ 6 h] < \ P {p - l)X E\v%\* + ip(p - l) P E\u?f H - 2 

Using that \u\ff < \u\ & H + t^j, we obtain by easy computations that 

j t EK\% +p V C x E\u!f H < \p(p-l){\ + e)E\v%\v + C{s,p, P ,C x ,v) 

for any e > (with the latter constant C being a suitable positive constant). If 
2vC x > Xq, we conclude as before by Gronwall Lemma that there exists p > 2 
such that 

E\ul l \ p H <C 7 Vt>0,Vn>l 
From now on, the proof proceeds as in the previous case. ■ 

4.1 Scaling 

Let us start with an heuristic digression. We recall that in the Kolmogorov 1941 
theory of turbulence (see, e.g., Sect. 6.3.1 in [10] dealing with the deterministic 
equations and [llj dealing with the stochastic equations), one believes that the 
following equality in law is approximatively true 

u(r + Xx) — u{r) m = w A 1 ^ 3 [u(r + x) — u{r)\ 

for any r, x S K 3 and for A in some range of small positive real numbers. 

(In the whole section, u(x), without the time variable, denotes a stationary 

solution.) 

This implies A" 1 / 3 [u(Xx) - u(0)] in = w u(x) - u(0). 

According to this result, we are interested in the scaled velocity u\, defined by 
the following scaling transformation 

u x {t,x) := A~ 1/3 u(A 2/3 t,A:r) 

for A G (0, 1); hence the function u\(t,x) is defined for x £ [0, ¥] . 

We assume that u = u(t,x) solves in the torus [0,L] 3 the modified Navier- 
Stokes equation (with $>(u) — ^|u| 4 u), with additive noise 

du + [-A$(u) + (u ■ V)it + X7q] dt = ^ a k ,jd0 ktj (t)hk,j 

(kj)eA 

(with respect to (fTOj) . there are the coefficients (Ju,j \ some of them may vanish 
and therefore we denote by A the set for the summation on cr^.j ^ 0. Condition 
© is satisfied if J2{k,j)eA \ a k, 3 \ 2 < oo.) 

The scaled velocity satisfies an equation very similar to this one. 
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Proposition 11 We have 

dux + [-A<f>(u x ) + (ux-V)ux + Vq\]dt= vuWkAWkj 

(fej)eA 

where q\ is a suitable function, h^ -{x) — hkj(Xx) and the processes 

fa (t) := X-y^ kJ (x^t) 

are independent standard Brownian motions. 

Proof. The rigorous proof has to be performed at the level of the integral 
weak formulation of the equation and it is tedious and elementary. We just 
point out the main (somewhat heuristic) arguments behind it. We have 



du x 
dt 



(t,x) 



dt 



( X 2 / 3 t,Xx) 



(u x • V) u x \ {t . x) = X 1/3 (u • V) u|( A2/3t>Ax ) 



dt 



(t) 



^.(t + dt)-^(t) 
dt 

>rl/3 ttA x2/3t + x2/3 *)-ttA x2/3t ) 



dt 



A l/3 
A l/3 



(3 kJ {X 2 'H + X 2 'Ht) - (3 kJ (X 2 / 3 t) 



X 2 / 3 dt 



dt 



(X 2 / 3 t) 



A^K)| (M) = AV3 A$(«)| (AV , tiAx) 



because 

D l \ux\ 4 ux = X-^ 3 D l 



X~ 5 / 3, i 



+ a 



-5/3 



(a 2 / 3 ;, 

a (x 2 '\ 



Ax 2 



= A • A" 5 / 3 



(A 2 / 3 t,Ax)| 2 ' 2 «(A 2 / 3 i,Ax) 

(A 2 / 3 t, Az) 2 u (A 2 / 3 ^ Ax) • A u (A 2 / 3 i, Ax) 
(A 2 / 3 t,Ax)" 

(X 2 / 3 t,\x) 



Xx 



D 



u2\u\ u- (D t u) + \u\ z ' z (Diu) 



X • A" 5 / 3 Di 



\u\ 4 u 



(\ 2 / 3 t,Xx) 
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and then 

A\u x fu x =A 2 -A^ 5/3 A[|u| 4 u 
(t,x) L 



(\ 2 / 3 t,\x) 



Now, for any r £ R 3 , define the (space) translation operator f as (fV)(x) = 
V(x + r), to be understood as an identity in the distributional sense. We say 
that a process V is spatially homogeneous if all the space increments SV(x, h) — 
V{x+h)— V(x) are statistically invariant with respect to the translation operator 

f: V(x + h)- V(x) ia = w V(x + h + r) - V{x + r). 

In the same way, we say that a process V is isotropic if the law of all the space 
increments 5V(x, h) do not change under simultaneous rotation 9 of the space 
variables and of the vector V. Since the space variable lives in a torus, only 
rotations of R 3 which leave the torus invariant are allowed. 

It is easy to check when the Wiener process on the r.h.s. of our equation 
enjoys these statistical invariances. Indeed 

E\w(t,x + h)-w{t,x)\ 2 =t Wk,j\ 2 \h k ,j(x + h)-h k ,j(x)\ 2 

(kj)eA 

If for any k with (k, j) £ A, each coefficient a k .j in front of cos(^fc • x) is equal, 
in absolute value, to a coefficient in front of sin(^-fc-x), then the second moment 

of the space increment is equal to tjp \°~k,j\ 2 \e l ~ k ' h — 1| 2 and therefore 

(fej)eA 

depends only on h. This implies that w is spatially homogeneous (considering 
the space variables in R 3 /[0, L] 3 ). On the other hand, 

E\0w(t,6x + 6h)-6w(t,6x)\ 2 =t ^ Wk^\ 2 \h kJ (9x + 9h) - h k ^(9x)\ 2 

(k,j)eA 

But h k j (9x + Oh) — h k ,j (9x) = hg-i k j(x+h) — h g -i k j(x). Then we can consider 
only rotations 9 such that (k,j) £ A (9k, j) £ A and in these cases w is 

isotropic if Icrgfcjl = Wk,j\ for all (k,j) £ A. 

Corollary 12 Let A be such that the process Y2(k j)eA a kd0k,j{t)hk t j(x) is spa- 
tially homogeneous and isotropic. For $(w) = ^|m| 4 m 7 consider the equation 

du + [-A$(u) + (u ■ V)u + \7q] dt = Y a k,jd/3h,j(t)hk,j(x) 

(fcj)GA (24) 



x £ [0, L 



3 



with initial velocity spatially homogeneous and isotropic, satisfying the assump- 
tions of Theorem \10\ and X)(fcj)eA l* 7 *^'! < 00 ■ Then there exists a solution u 
spatially homogeneous and isotropic for any t > 0. For any %j) £ T>°° , we have 

E[\(u(t, Xe) - u(t, 0), V) | p ] = \ p/3 E[\ (u x (\- 2 'H, e) - u x (\- 2 'H, 0), $ |*] (25) 
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where u\(t,x) is spatially homogeneous and isotropic, and solves the equation 

du x + [-A^(ux) + (ux-V)u x +Vqx]dt = ^ a kJ d0^ d (t)h k j(>a) 

(fej)eA 

xe[0,±] 3 

(26) 

with initial velocity u\(0, •) = A 1 / 3 u(0, A-). 

Remark 13 The statistical invariance for the solution is obtained from the 
same property of the Galerkin approximations, as done in a similar context in 
U8\j . Indeed, we construct a solution as limit of a Galerkin subsequence. But 
the statistical invariance for the Galerkin processes u n (for any n) is easy to 
show, since for any n the finite- dimensional problem has a unique solution. 
Notice that in our case we can trivially consider a vanishing initial velocity. 

Now, let us consider the structure function of order p (p = 1,2,...) with 
respect to a stationary solution u of the modified Navier-Stokes equation ([24)) : 

S p {\) :=E[\u(Xe)-u(Q)\ p } 

(e is a unitary vector in R 3 and A <E (0, 1)). Similarly, we can work with the 
longitudinal structure function. 

We point out that Theorem [TO] provides the existence of a stationary solution 
leaving in X, but this is not enough to define the velocity in every point of the 
torus. We would need to analyze the regularity of stationary solutions, but we 
decide to postpone the study of existence of more regular stationary solutions 
to future work (it would be enough to have the law of u supported by the space 
C°(T)). 

According to the previous Corollary, for the structure function we get that 

S p (\) = X p / 3 E[\u x (e)~ux(0)\ p } (27) 

Kolmogorov 1941 theory states (see, e.g., [TO], Sect. 6.3.1) that 

S p (X) = C p e p/3 X p / 3 (28) 

where C p are dimensionless and e is the mean energy dissipation rate. For 
p = 2, l|28p is the so called two-thirds law of turbulence, which is supported by 
experimental results. For p — 3, ((28)) is the four-fifths law of turbulence (C3 = 
— |), deduced from the assumptions of homogeneity, isotropy and finiteness of 
the energy dissipation. For p > 3, (|28[) is not confirmed by experimental data 
and its truthfulness is questionable. 

According to the above Proposition and Corollary, we shall provide a rela- 
tionship similar to ([28)) for our model ([24)) . 

Keeping in mind ([27)) . we investigate the behaviour of E ^ux(e) — ma(0)| p ] in 
order to get insights on the structure function. First, we remark that (3k ,j (t) and 
Pkj W are unitary Brownian motions. Thus, the random forces in equations 
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(pM)) and (|26[) are the same in law; what changes is the dimension of the torus and 
correspondingly the eigenvectors hkj. No other terms in the equation (when 
projected onto H) change with the scaling. It is important to point out that 
this property is not true for the usual Navier-Stokes equation or for the Prouse 
model introduced at the beginning; namely, after the scaling the viscous term 
vAu becomes \~ 4 ^ 3 i>Au\ (and then the two problems for u and u\ are very 
different, because the scaled viscosity v\ = \~^l z v explodes as A — > 0). 

Since the coefficients a^j in the noise do not change with the scaling transfor- 
mation, the mean energy introduced in a unit of volume per unit of time by the 
stochastic forcing term is independent of A and is equal to j*g J2(k j)eA l CT *:.j| 2 - 
When A — > 0, the size of the domain becomes bigger and bigger but the unitary 
energy does not change. 

If, as in turbulence theory, we assume that during the motion there is energy 
transfer from large scales to small scales with a universal cascade mechanism 
depending only on the unit volume energy, then we would conclude that any 
stationary state is independent of A. Hence E \\u\(e) — u\(fS)\ p \ = k p for any A. 
Coming back to (|2"5"|) . we conclude that 

s p (\) = k p y/ z 

for any < A < 1 and p eN. 

Summing up, we have the following result, providing a result on the struc- 
ture function (of any order p) under two assumptions. The first assumption is 
technical and can be removed as soon as we are able to prove existence of reg- 
ular stationary solutions. On the other hand, the second assumption on energy 
cascade has to be considered as an hypothesis quite hard to justify rigorously 
(as it is for fluids modeled by the Navier-Stokes equations). 

Claim 14 Let us assume that system (|24p has a stationary solution u, which 
at any fixed time has law supported by the space C'° (T) . 

Let us further assume that there is energy transfer from large scales to small 
scales with a universal cascade mechanism depending only on the unit volume 
energy. 

Then for the structure function, given any < A < 1 and p e N we have 

S P (X) = k p X p / 3 
for some constant k p independent of A. 

5 Appendix 1: a priori estimates 

We present the estimates on the Galcrkin approximations; these are quite stan- 
dard (see, e.g., [9] and [H]). Besides the usual estimates we need 
also (|37j) to prove uniqueness. 

Let (it")t>o be a continuous adapted solution of equation (Tl3|) . Let 

r£ = inf{t>0: \u\% = R] 
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We have 



//,: + / - tt„B«, uj)] ds + / 7r„G«)dw(s) 

o ,/o 



It, 







-A$« AT „) - 7r nJ B« Ar „,< Ar „) 



l{ S <r^}rfs 



For p > 2 apply Ito formula to |m™ At ,i|^: 

d\ U ? AT n\ P H <p|< AT g|^ 2 « T3)( i< AT g)H 



+ |p(p- 1)K T »|S 2 |k„G«^)llH5(H)l{ t <^}^ 



Then, integrating in time, we have 



,n ip 



nip 



ft 

,n IP— 2 



+ P / KatjIh (-^«Arg)-7rn S «ATg 5 <Ars)><Ar«)H 1 { S <r S }rfs 



+ P / Kat«Ih 2 «ArS' 1 {^<rS} 7r ™G'« ATS )dw(s))_ff 



1 



+ 2^-!) / l«?Argl£" ^G« AT „) 
Then, by Lemma [1] and (@| 



HS(H) 



1 {s<T%}ds 



n ip— 2|| n 



i (is 



< / l< ArS l?T 2 « ATS ,i {s <r S }^GK AT „)^( s )>H 

+ |p(P" 1) f l<Ar S llr 2 ||G« Ar g)||^ (H) l {s < T g } d S 



i) 



By assumption © 



|U?A T -ISr+^ / l«?Ar S iri«?Ar 5 llvl{.<T S }<k 

^>(P - 1) j[* Ka. S I?T 2 (Ao|< AtS | 2 h + p) l { .<r S }da 

(29) 



IP-2, 



<KI?r+P 



where 



M?= / |< Ar ,|?T 2 (< Ar ,,l {s < TS } 7 r„G(< Ar „)^( S )) ff 

M 
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is a square integrable martingale. 
Therefore 



sup K ATg |^ < K\ p H +p sup 
t£[0,r] te[o,r] 

1 



+ ^P(p - !)(Ao + P) / l<AT«lif 1 {s<Tg}rfs + _ l )P r 

(30) 

By Burkholder-Davis-Gundy inequality, we estimate the supremum of the mar- 
tingale M t ra ; for some constant C > we have 



pE sup 

o<to 



< CpE ^ |< ArS |^ 2 ||G« Ar „)||^ (ff) l {s<rS} d S ) 



1/2 



Then by assumption 



pE sup 

te[o,r] 



Ml 



< E 



sup K T j^ a cpi/ |tt? AT »l5r a (Ao|u; Ar5 |^+p)i {i<T » } d» 

0<t<r . / . . 



1/2' 



< ( E sup 

0<t< 



J< Ar g|^ ( EC2P2 [ K^ R \ P H\^K^n\ 2 H + P )l {s<r n } ds\ 

< l -E sup K rg |£ + icV(Ao+p)E ri<Ar S l?rl{ s <r S }rfs+icVpr 

^ 0<t<r R Z J q A 



(31) 



Then, by (30]) and (J31 
' - 1 ,n ip 



£ sup |u? AT »|» <SK|^ + -CV(A + P) / S|< AT „|? r l {5<T » } d s + -CVpr 



te[o,r] 



<SK|^ + icV(A +p) f 25 sup |< |^ S + ic 2 P V 
* Jo te[0,a] ^ 

(32) 



By Gronwall lemma, for any r > we have 

E sup K ATJ ,|£ <d 

0<i<r K 



(33) 



for some positive constant C\ — Ci(p,T, Xq, p,m p ) independent of n and R. 
Here m p = E\uq\ p h . Notice that E [\uq \ p h ] < m p . 



27 



Coming back to (|29[) . with similar arguments we also obtain 

1 s/\t r \'h \\ u s/\Trt \\v ^{s<t r } 

for a new positive constant C depending on m p ,p, Ao, p, T but not on n, R. For 
p = 2we have 

e[ \K ATa \\ 2 v l {s<TR} ds < C 2 (34) 
Jo 

for some positive constant C 2 = C 2 (T, Ao, p, 7712) independent of n and R. 
For the last estimate, we proceed as follows. First, from @ we have that 

u(x)\ 1+b dx = \u(x)\ 1+b l { \ u ( x )\< K} dx + / |w(x)| 1+ ''l{| u ( x )| > _ R -}dx 
r jt jt 

<if 1+b |T| + / \u(x)\ 2 — a(\u(x)\)dx 
Jt a i 

= K 1+b \T\ + -{<5>{u) 1 u) H 

i.e. 

(^W.^^axK+i-ai^lTI (35) 

Moreover 

2\{B{u,u),A- 1 u)\ < 2|tt|£*|«| V / < H^ + H^<ai|u|^+ 6 & + C' + C , |u|l f (36) 
Apply Ito formula to \u r t L AT J\ 2 v , = (< AT g, A^u^^h and get 

r-t 



r> 



|< AT „|^ =Klv' -2 / ($« Ar „),< Ar „)Hl {s <. s} d S 
c 

(Tn-B(u7 ATg ,ttJ VTg ),i4- 1 «; AT5 ) H l {s<T « } ds 
+ m + f \WnA~ 1/2 G(u^ AT n)\\ 2 HS{H) l {s<T n } ds 



where 



Aft = 2 



«ArS- l{ S <r«}7Tn^ 1 G« AT „)d U ;( S ))H 



is a square integrable martingale. 
Use ([35J and (011); then 



E\u^ Ar n\ 2 vl + ai Ej 



T 

U sAt£ \~L 1 + b ^{s<r™}d>S 



<E\uX>+CTE sup \u? AtR \ 2 h 

0<t<T H 



CXoE^ \u: AT , A \ 2 H l {s<T . } ds + C(l + P )T 
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According to (|3"3"|) (for p = 2), we conclude that 

E f T \K AT ^ b l {s <^}ds<C 3 (37) 
for some positive constant C3 = Cs(T, Ao, ai, Ci, C2). 

6 Appendix 2: a compactness result 

Let X be the closure of 2?°° w.r.t. the norm 

3 



1/6 



u| x := / {\u{x)\ A \Vu(x)\ 2 +4\u(x)\ 2 J2Mx) ■ d iU (x)] 2 }dx 



r i=i 



Theorem 15 X C L 6 n H and the immersion is compact. 
Proof. First observe that for smooth fields u we have 
di {\u\ 2 u) = \u\ 2 diu + 2u(u ■ diu) 

Hence 

\\\u\ 2 uf v <cj \u\ 4 \Vu\ 2 dx < £7|«|| 
and thus by Poincare inequality 

|«|£. = |N 2 «£,<||H 3 U ||* <C|«|| 

This proves that the closure of T>°° with respect to the L 6 -norm is a space bigger 
than its closure with respect to the X-norm; hence X C L 6 n H . 

Moreover, if {u n } is a bounded sequence in X, it is bounded in L 6 and 



{ |u n | 2 it„ } is also bounded. By Rellich Theorem, the sequence {|u„| z u„} 

is relatively compact in L 2 and so there exists a subsequence {\u nk \ 2 u„J con- 
verging strongly in L 2 to some field £; we also have that {u nk } converges weakly 
in L 6 to some field u. The strong convergence implies in particular that 



2 2 



U 



I 6 

n k lL 6 



\ u n k \ u n k 



Thus, if we prove that |£|£ 2 = |u|^ 6 , then from the weak convergence of {u nk } 
to u in L 6 and the convergence of norms |iinJ]L6 — > M L6 , we deduce that \u nk ] 
converges strongly to u in L 6 and the proof of the compact embedding will be 
complete. 

So it remains to show that |£|t 2 = \uUe- Let us introduce the function 

,(*) = { ^ if . 

I if C (x) = 
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Let us prove that there is a subsequence {u n > k } such that 

u, n i — > v a.s. on T. 



(38) 



This implies v — u (the a.s. limit and the L 6 weak limit must coincide, since 
by Vitali theorem there is strong convergence in any L p with p < 6). Since 
\v (x)\ e — |£ (x)| 2 where £ (a;) ^ 0, we have \£,\^2 — Ml 6 — M^e, as we want. 

Thus it remains to prove (|3"8|) . The strong convergence above implies that 
there is a subsequence {|u„'J 2 u„j_} that converges to £ a.s. on T. Let x 6 T 
be such that |u„' (x)| 2 w n ' (x) — * £(a;). Taking the norm in R 3 , this implies 
that \u n i (x)\ 3 — > |£(a:)|, hence |u n / (x)| — > ^(x)! 1 / 3 . If £(x) = 0, this implies 



u n' k (x) —> 0, as we want in ([38)) . If £ (x) =^ 0, this implies 
and 



(;,) 



7^ eventually 



(x) 



V (») u n' (x) 



(X) 



\H X )\ 



2/3 ■ 



Thus (|3"5|) is true. The proof is complete. 
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